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Abstract 

We calculate the branching ratios and CP asymmetries for B° — > vr + p~, B° — ► p + vr _ , 

^ . B + -> p + 7T°, 5+ -> vr+p , 5° -> vr°p , -> tt+cj and 5° -> vr°u; decays, in a pertur- 

H ' 

bative QCD approach. In this approach, we calculate non-factorizable and annihilation 
type contributions, in addition to the usual factorizable contributions. Our result is in 
agreement with the measured branching ratio of B°/B° — > ^p T , B ± -> n ± p°, n ± u by 
CLEO and BABAR collaboration. We also predict large CP asymmetries in these decays. 
These channels are useful to determine the CKM angle <fo- 
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1 Introduction 



The rare decays of B mesons are getting more and more interesting, since they are useful for 
search of CP violation and sensitive to new physics. The recent measurement of B — > Tip 
and ttuj decays by CLEO Collaboration M arouse more discussions on these decays 0. The 
B — > up, huj decays which are helpful to the determination of Cabbibo-Kobayashi-Maskawa 
(CKM) unitarity triangle 02 have been studied in the factorization approach in detail || [|. 
In this paper, we would like to study the B — > irp and txuo decays in the perturbative QCD 
approach (PQCD), where we can calculate the non-factorizable contributions as corrections to 
the usual factorization approach. 

In the B — > 7rp, txlo decays, the B meson is heavy, sitting at rest. It decays into two light 
mesons with large momenta. Therefore the light mesons are moving very fast in the rest frame 
of B meson. In this case, the short distance hard process dominates the decay amplitude. The 
reasons can be ordered as: first, because there are not many resonances near the energy region 
of B mass, so it is reasonable to assume that final state interaction is not important in two-body 
B decays. Second, With the final light mesons moving very fast, there must be a hard gluon 
to kick the light spectator quark (almost at rest) in the B meson to form a fast moving pion 
or light vector meson. So the dominant diagram in this theoretical picture is that one hard 
gluon from the spectator quark connecting with the other quarks in the four quark operator of 
the weak interaction. There are also soft (soft and collinear) gluon exchanges between quarks. 
Summing over those leading soft contributions gives a Sudakov form factor, which suppresses 
the soft contribution to be dominant. Therefore, it makes the PQCD reliable in calculating 
the non-leptonic decays. With the Sudakov resummation, we can include the leading double 
logarithms for all loop diagrams, in association with the soft contribution. Unlike the usual 
factorization approach, the hard part of the PQCD approach consists of six quarks rather than 
four. We thus call it six-quark operators or six-quark effective theory. Applying the six-quark 
effective theory to B meson decays, we need meson wave functions for the hadronization of 
quarks into mesons. All the collinear dynamics are included in the meson wave functions. 

In this paper, we calculate the B — > it and B — > p form factors, which are input parameters 
used in factorization approach. The form factor calculations can give severe restrictions to 
the input meson wave functions. We also calculate the non-factorizable contributions and the 
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annihilation type diagrams, which are difficult to calculate in the factorization approach. We 
found that this type of diagrams give dominant contributions to strong phases. The strong 
phase in this approach can also be calculated directly, without ambiguity. In the next section, 
we will briefly introduce our method of PQCD. In section 3, we perform the perturbative 
calculations for all the channels. And we give the numerical results and discussions in section 
4. Finally section 5 is a short summary. 



2 The Frame Work 

The three scale PQCD factorization theorem has been developed for non-leptonic heavy meson 
decays H. The factorization formula is given by the typical expression, 



C(t) x H(x,t) x x exp 



-s(P,b)-2 f ^l q (a s m 
Jl/b u 



(1) 



1/6 fJ, 

where C(t) are the corresponding Wilson coefficients, are the meson wave functions. And 
the quark anomalous dimension 7 g = —a s /ir describes the evolution from scale t to 1/6. 

Non-leptonic heavy meson decays involve three scales: the W boson mass mw, at which 
the matching condition of the effective Hamiltonian are defined, the typical scale t of a hard 
subamplitude, which reflects the dynamics of heavy quark decays, and the factorization scale 
1/6, with b the conjugate variable of parton transverse momenta. The dynamics below 1/6 
scale is regarded as being completely non-perturbative, and can be parameterized into meson 
wave functions. Above the scale 1/6, PQCD is reliable and radiative corrections produce two 
types of large logarithms: ln(raw/£) and ln(t6). The former are summed by renormalization 
group equations to give the leading logarithm evolution from mw to t scale contained in the 
Wilson coefficients C(t). While the latter are summed to give the evolution from t scale down 
to 1/6, shown as the last factor in eq.(jl]). 

There exist also double logarithms ln 2 (P6) from the overlap of collinear and soft divergence, 
P being the dominant light-cone component of a meson momentum. The resummation of 
these double logarithms leads to a Sudakov form factor exp[— s(P, b)], which suppresses the 
long distance contributions in the large 6 region, and vanishes as 6 > 1/Aqcd- This factor 
improves the applicability of PQCD. For the detailed derivation of the Sudakov form factors, 
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see ref.p, [7|. Since all logarithm corrections have been summed by renormalization group 
equations, the above factorization formula does not depend on the renormalization scale \x. 

With all the large logarithms resummed, the remaining finite contributions are absorbed 
into a hard sub-amplitude H(x, t). The H(x, t) are calculated perturbatively involving the four 
quark operators together with the spectator quark, connected by a hard gluon. When the end- 
point region (x — > 0, 1) of wave function is important for the hard amplitude, the corresponding 
large double logarithms a s ln 2 x shall appear in the hard amplitude H(x,t), which should be 
resummed to give a jet function St(x). This technique is the so-called threshold resummation ||. 
The threshold resummation form factor S t (x) vanishes as x — > 0, 1, which effectively suppresses 
the end-point behavior of the hard amplitude. This suppression will become important when 
the meson wave function remains constant at the end-point region. For example, the twist-3 
wave function <fi^ and (fa are such kinds of wave functions, which can be found in the numerical 
section of this paper. The typical scale t in the hard sub-amplitude is around y/KMs- It is 
chosen as the maximum value of those scales appeared in the six quark action. This is to 
diminish the a 2 corrections to the six quark amplitude. The expression of scale t in different 
sub-amplitude will be derived in the next section and the formula is showed in the appendix. 



2.1 Wilson Coefficients 

First we begin with the weak effective Hamiltonian H e ff for the AB = 1 transitions as 

10 

v ub v: d {c x oi + c 2 od - v tb v; d £ a o t 



Keff ~ — 



V2 



i=3 



(2) 



We specify below the operators in H. e ff for b — > d: 



0\ 


= dofl^Lup ■ 


upjpLba , 




= daJ^LUa ■ 


up'Jt+Lbp , 




= da^Lba ■ 




o 4 


= daTpLbp ■ 


Eg' q'f3l»Lq' a , 




= doTf^Lba ■ 




o 6 


= dcl^Lbp ■ 


Eg' q'pl^Wa > 


7 


= \d a ^Lb a 


■ Eg' e q 'qpi»Rqp , 


o 8 


= \d a ^Lb p 


■ Eg' eg'q'p^Rq'a 


o 9 


= \d a ^Lb a 


■ Eg' egiq'pl^Lq'p , 


010 


= \d a ^Lb p 


■ Eg' z q >q'f3i»Lq' a 



(3) 



Here a and (3 are the SU(3) color indices; L and R are the left- and right-handed projection 
operators with L — (1 — 75), R — (1 + 75). The sum over q 1 runs over the quark fields that are 
active at the scale \i = 0(m&), i.e., (q'e{u, d, s, c, b}). 



The PQCD approach works well for the leading twist approximation and leading double 
logarithm summation. For the Wilson coefficients, we will also use the leading logarithm 
summation for the QCD corrections, although the next-to-leading order calculations already 
exist in the literature || . This is the consistent way to cancel the explicit \x dependence in the 
theoretical formulae. 

If the scale nib < t < mw, then we evaluate the Wilson coefficients at t scale using leading 
logarithm running equations |J, in the appendix B of ref. [10]. In numerical calculations, we 
use a s = Ait /[f3\\n(t 2 / K^q CD )] which is the leading order expression with A { q ] cd = 193MeV, 
derived from ^q CD = 250MeV. Here /3\ = (33 — 2n/)/12, with the appropriate number of active 
quarks rif. rif — 5 when scale t is larger than m&. 

If the scale t < nib, then we evaluate the Wilson coefficients at t scale using the formulae in 



appendix C of ref. [10| for four active quarks (rif — 4) (again in leading logarithm approxima- 
tion). 



2.2 Wave Functions 



In the resummation procedures, the B meson is treated as a heavy-light system. In general, 
the B meson light-cone matrix element can be decomposited as |TI], [TI| 



iki z 



o (2tt) 4 
i I 



■ z (0\b a (0)d p (z)\B(p B )) 



V2N C 



0>B + m B )-f 5 



(4) 



Pa 



where n = (1,0, Ot), and v = (0,1, Ot) are the unit vectors pointing to the plus and minus 
directions, respectively. From the above equation, one can see that there are two Lorentz 
structures in the B meson distribution amplitudes. They obey to the following normalization 
conditions 



d 4 h 



(2tt) 4 ^ Bv " 1/ 2v/2iV7 J (2tt) 4 
In general, one should consider both these two Lorentz structures in calculations of B meson 



d 4 h 



>B 



0. 



(5) 



decays. However, it can be argued that the contribution of 0^ is numerically small |I3" |, thus 
its contribution can be numerically neglected. Therefore, we only consider the contribution of 
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Lorentz structure 



B 



V2Nc 



Jb + mB)750s(ki^ 



(6) 



in our calculation. We keep the same input as other calculations in this direction 10, 13, 14 



and it is also easier for comparing with other approaches []T2], |i5| . Through out this paper, we 
use the light-cone coordinates to write the four momentum as (kf, k^ , k±). In the next section, 
we will see that the hard part is always independent of one of the kf and/or /cf , if we make 
some approximations. The B meson wave function is then the function of variable (or kf) 
and k±. 

MK, *t) = I dkf^kt, fcr, ki). (7) 



The 7r meson is treated as a light-light system. At the B meson rest frame, pion is moving 
very fast, one of kf or ki is zero depends on the definition of the z axis. We consider a pion 



moving in the minus direction in this paper. The pion distribution amplitude is defined by [1C 



<tt-(p)KMMo)|o> 



V2N c Jo 



7 5 jP(j) n (x) + m O 7 5 0p(x) - moa^^P^ 
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(8) 



For the first and second terms in the above equation, we can easily get the projector of the 
distribution amplitude in the momentum space. However, for the third term we should make 
some effort to transfer it into the momentum space. By using integration by parts for the third 
term, after a few steps, eq.(H) can be finally changed to be 



<n-{P)\d a {z)u p (0)\0> 



ixP-z 



75 



(x) + m O 7 5 0p(x) + m [7 5 (/; jh - 



/3a 



(9) 



where = and vector v is parallel to the ir meson momentum p n . And m = 

m^j(m u + raj) is a scale characterized by the Chiral perturbation theory. In B — > up decays, 
p meson is only longitudinally polarized. We only consider its wave function in longitudinal 
polarization [|13|, [17]]. 

1 rl 



<p-(P,e L )\d a (z)up(0)\0>-- 



{4 [Mfa) + ™ P <P P (x)} +m p <f) s p {x)} (10) 



y/2N~ c J0 

The second term in the above equation is the leading twist wave function (twist-2), while the 
first and third terms are sub-leading twist (twist-3) wave functions. 
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The transverse momentum k L is usually conveniently converted to the b parameter by 
Fourier transformation. The initial conditions of 4>i(x), % — B, tt, are of non-perturbative 
origin, satisfying the normalization 

j^ l ( x ,b = 0)dx = ^=f l , (11) 

with fi the meson decay constants. 

3 Perturbative Calculations 

In the previous section we have discussed the wave functions and Wilson coefficients of the 
factorization formula in eq.(|l]). In this section, we will calculate the hard part H(t). This 
part involves the four quark operators and the necessary hard gluon connecting the four quark 
operator and the spectator quark. Since the final results are expressed as integrations of the 
distribution function variables, we will show the whole amplitude for each diagram including 
wave functions. 

Similar to the B — > hit decays [fTO , there are 8 type diagrams contributing to the B — > up 



decays, which are shown in Figure 1. Let's first calculate the usual factorizable diagrams (a) 
and (b). Operators Oi, 2 , 3 , 4 , 9 , and O w are (V — A)(V — A) currents, the sum of their 
amplitudes is given as 

F e = sV^irCFGFfpfTipm^e-p^) / dx\dx 2 / bidbib 2 db 2 (f)B(xi,bi) 

Jo Jo 

{ [(1 + x 2 )<t>i(x 2 , b 2 ) + rv(l - 2x 2 ) (x 2 , b 2 ) + C(x 2 , 6 2 ))] a s {t\) 
h e (x!,x 2 , h, b 2 ) exp[-S ab (tl)} 

+2r^(x 2 ,b 2 )a s (t 2 e )h e (x 2 ,x 1 ,b 2 ,b 1 )exp[-S ab (t 2 e )]} , (12) 

where r n = m /m B = m^l\mB{ra u + m d )}; Cp = 4/3 is a color factor. The function h e , the 
scales f e and the Sudakov factors S ab are displayed in the appendix. In the above equation, 
we do not include the Wilson coefficients of the corresponding operators, which are process 
dependent. They will be shown later in this section for different decay channels. The diagrams 
Fig.l(a)(b) are also the diagrams for the B — >• it form factor Ff^ n . Therefore we can extract 
pB^ir from eq.flTJ). 

Fftf = 0) = / e -. (13) 

V2& F f p m p (e-p 7T ) 
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The operators O5, Oq, O7, and Os have a structure of (V — A)(V + A). In some decay channels, 
some of these operators contribute to the decay amplitude in a factorizable way. Since only the 
vector part of the (V+A) current contribute to the vector meson production, 



{tt\V - A\B)(p\V + A\0) = (ir\V - A\B)(p\V - A\0), 



(14) 



the result of these operators is the same as eg . (|i~2"D . In some other cases, we need to do Fierz 
transformation for these operators to get right color structure for factorization to work. In 
this case, we get (S-P)(S+P) operators from (V-A)(V+A) ones. Because neither scaler nor 
pseudo-sealer density gives contribution to a vector meson production (p\S + P\0) = 0, we get 



F 



0. 



(15) 



6 


P 

V 
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(a) 





Figure 1: Diagrams contributing to the B — > irp decays (diagram (a) and (b) contribute to the 
B — > it form factor). 

For the non-factorizable diagrams (c) and (d), all three meson wave functions are involved. 
The integration of 63 can be performed easily using 5 function 5(63 — 61), leaving only integration 
of b\ and 62- For the (V — A)(V — A) operators the result is 



32 

y 



•f'2 



V3irCFGFTn p m^(e ■ p n ) / dxidx 2 dx 3 / bidbib 2 db2 4>B(xi,bi) 

Jo Jo 

4>„(x 2 , h) - 2r n (j)°(x 2 , h) 4> p (x 3 , b 2 )h d (x 1 ,x 2 , x 3 , b ly b 2 ) expf-S^^)] . (16) 



S 



For the (V — A)(V + A) operators the formula is different, 

64 i— r°° 
M e = — V3irCFGFrn 2 m B (€ ■ p n ) / dx\dx 2 dx 3 / bidbib 2 db 2 4>b{ x 1i &i) x 
3 jo jo 

{r w (s 3 - x 2 ) [<f>%(x 2 , &i)0*(x 3 , 62) + C(^2, bx)(j) s p {x^ 6 2 ) 

-7V(x 2 + ^3) [^(»2, &l)0pO&3, 62) + CO^, &l)0p(x 3 , 62) 

+x 3 0^(x 2 , 61) P (x 3 , 6 2 ) - p (x 3 , 6 2 ) } h d (x u x 2 , x 3 , 61, 6 2 ) expf-'S'ajfo)] . (17) 

Comparing with the expression of M e in eq.([L6]), the (V— A)(V+A) type result is suppressed 
by m p /m B . 

For the non-factorizable annihilation diagrams (e) and (f), again all three wave functions 
are involved. The integration of b 3 can be performed easily using S function S(b 3 — 62)- Here 
we have two kind of contributions, which are different. M a is contribution containing operator 
type (V — A)(V — A), and is contribution containing operator type (V — A)(V + A). 
32 1- r 1 r°° 

— v3vrCi?Gi?m p m^(e ■ p w ) / dx\dx 2 dx 3 / b\db\b 2 db 2 4>b( x i, 61) x 
3 Jo jo 

{ x 2 (f)^(x 2 , b 2 )<j) p (x 3 , b 2 ) + r % r p (x 2 - x 3 ) (^(x 2 , & 2 )0 P (x 3 , b 2 ) + C(x 2 , 6 2 )0p(x 3 , 6 2 )) 
+r ff r p (a;2 + x 3 ) fe(x 2 , & 2 )</> P (a: 3 , 62) + ^(^2, b 2 )(f) s (x 3 , b 2 



Mi 



xh 1 f (x 1 ,x 2 , x 3 , 61, 6 2 ) exp[-5 e/ (t})] 

- x 3 <f>^(x 2 , 62)^(3*, 62) + r w r p (x 3 - x 2 ) (^(x 2 , fr 2 )0 P (x 3 , b 2 ) + C(x 2 , b 2 )(f) s p (x 3 , b 2 
+r 7r r p (2 + x 2 + x 3 )4>^(x 2 , b 2 )4> p (x 3 , b 2 ) - r n r p (2 -x 2 - x 3 )(f)°(x 2 , b 2 )(f) p (x 3 , b 2 ) 
h? f (xi,x 2 , x 3 , 61, 6 2 ) e w[-Sef{t})]} , 

— — v3TrCFGF'rn p m B (e-p 7T ) / dx\dx 2 dx 3 \ b\db\b 2 db 2 4>b ( x i, b\) x 



(18) 



{ [ararv^p (ar 3 , & 2 ) (<^(x 2 , 6 2 ) + 6 2 )) - x 3 r p <f4(x 2 , h) (f p {x 3 , b 2 ) + p (x 3 , b 2 

xhj(xi,x 2 , x 3 , 61, 6 2 ) exp[-5 e/ (t})] 

(2 - x 2 )r 7r p (x 3 , 6 2 ) (4>n( x 2, h) + K( x 2i b 2 



(19) 

-r p (2 -x 3 )0^(x 2 ,& 2 ) (0 P (x 3 ,6 2 ) + p (x 3 ,6 2 ))] /ij(xi,x 2 ,x 3 ,6i,6 2 ) expf-S^tJ)]} , 

where r p = m p /mg. The factorizable annihilation diagrams (g) and (h) involve only ir and p 
wave functions. There are also two kinds of decay amplitudes for these two diagrams. F a is for 
(V — A)(V — A) type operators, and is for (S — P)(S + P) type operators. 

F a = 8V2CFGF^fBT r n P 'mB( e 'PTr) / dx\dx 2 \ b\db^b<> 

Jo Jo 
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[ x 2 (j)^(xi, bi)(J) p (x 2 , b 2 ) - 2(1 - x 2 )r 7r r p (j)^(x 1: &i)</> P (x 2 , b 2 ) 

+2(1 + x 2 )r 7r r p 0^(xi, bi)4> p (x 2 , b 2 ) a s (t\)h a (x 2 , x iy b 2 , b^ exp[-S gh (tl)} 

- Xi(j)£(xi, b 1 )(p p (x 2 , b 2 ) + 2(1 + x 1 )r 7r r p (j)P(x 1 , b 1 )(p s p (x 2 , b 2 ) 

-2(1 - xiY^p^xi, bx)(j) s p (x2, b 2 ) a s (t 2 e )h a (xi, x 2 , b lt b 2 ) exp^S^^)]} 

16-v/2Ci?Gi?7r/Bm p m|(e • p^) / dx\dx 2 \ b\db\b 2 db 2 

Jo Jo 

x | 27v</>f(xi, b 1 )(j)p(x 2 , b 2 ) + x 2 r P 4>^(xi, h) (4> p (x 2l b 2 ) - 4> l p (x 2 , b 2 )^ 
a s (t l e )h a (x 2 , Xt, b 2 , 6i) exp[-S gh (t l e )} 



+ 



x x r n (^(xijfti) -0^(xi,6i)j 4> p (x 2 ,b 2 ) + 2r p (j) 7r (xi,b 1 )(j) s (x 2 ,b 2 



a s (tl)h a (x 1} x 2 , bi, b 2 ) exp[-S gh (t 2 e )]} , 



(20) 



(21) 



In the above equations, we have used the assumption that X\ << x 2 ,x$. Since the light 
quark momentum fraction X\ in B meson is peaked at the small region, while quark momentum 
fraction x 2 of pion is peaked around 0.5, this is not a bad approximation. The numerical results 
also show that this approximation makes very little difference in the final result. After using this 
approximation, all the diagrams are functions of k^ = of B meson only, independent 

of the variable of kf. Therefore the integration of eq.(^) is performed safely. 
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Figure 2: Diagrams contributing to the B — > irp decays (diagram (a) and (b) contribute to the 
B — > p form factor Aq^ p ). 



If we exchange the ir and p in Figure 1, the result will be different for some diagrams. 
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Because this will switch the dominant contribution from B — ► 7r form factor to B — > p form 
factor. The new diagrams are shown in Figure 2. Inserting (V — A)(V — A) operators, the 
corresponding amplitude for Figure 2(a) (b) is 



ep 



S\/2nC F G F f v m p m 2 B {e ■ p^) / dx\dx 2 / bidbib 2 db 2 <p B (xi,bi) 

Jo Jo 

x { [(1 + x 2 )0 P (x 2 , b 2 ) + (1 - 2x 2 )r p (f p {x 2 , b 2 ) + 0*(x 2 , b 
a s (t\)h e (xi, x 2 , bi, b 2 ) exp[-S ab (tl)} 



+2r p <f) s (x 2 , b 2 )a s (t 2 e )h e (x 2 , x ± , b 2 , h) exp[-S ab (t 



(22) 



These two diagrams are also responsible for the calculation of B — > p form factors. The form 
factor relative to the B — > 7rp decays is Aq^ p ', which can be extracted from eq. ([22]) 



= 0) 



ep 



V2G F f ir m p (€ ■ p n ) 
For (V — A)(V + A) operators, Figure 2(a) and 2(b) give 

Ff p = -16V27rC F G F f 7T m p r 7r m 2 B (e ■ p n ) dx\dx 2 b\dbib 2 db 2 (j)B{xi,bi) 
x { 4> P (x 2 , b 2 ) - r p x 2 4> p (x 2 , b 2 ) + (2 + x 2 )r p (f) p (x 2 , b 2 ) 
xa s (tl)h e (x u x 2 , h, b 2 ) exp[-S ab (tl)] 

+ Xi<j) p (x 2 ,b 2 ) + 2r p <p s p (x 2 ,b 2 ) a s (t 2 e )h e (x 2 , x x , b 2 , h) expf-S^t*)]} 

For the nonfactorizable diagrams Figure 2(c) and 2(d) the result is 

32 /— f 1 f°° 

M Pn = — —y3TrC F G F m p m 2 B (e-p 7r ) / dx\dx 2 dx 3 / bidbib 2 db 2 (pB(xi,bi) 

Jo Jo 



(23) 



(24) 



XX 2 



(p P (x 2 ,b 2 ) -2r p (p t Jx 2 ,b 2 ) (p„(x 3 , b^h^xt, x 2 , x 3 , b u b 2 ) exp[-S cd (t d )} . (25) 



For the nonfactorizable annihilation diagrams (e) and (f), we have M ap for (V — A)(V — A) 
operators and for (V — A)(V + A) operators. 

32 r- f 1 f°° r 

M ap = — V3nC F G F m p m 2 B (e-p 7T )j^dx 1 dx 2 dx 3 b x db x b 2 db 2 <j) B (xi, h) |exp[— 5 e /(t})] 

x 2 <f)^{x 3 , b 2 )4> p (x 2 , b 2 ) + r n r p (x 2 - x 3 ) (^(x 3 , & 2 )0*(x 2 , 62) + b 2 )<f) p {x 2 , b 2 )) 

+r 1T r p {x 2 + z 3 ) (0^X3, b 2 )(f) p (x 2l b 2 ) + 0^(x 3 , b 2 )(j) s p (x 2l 6 2 )) h/fai, x 2 , z 3 , &i, 62) 

x 3 ^(a;3, 6 2 )0 p (x 2 , 6 2 ) + r n r p (x 3 - x 2 ) U%(x 3 , b 2 )(j) P (x 2 , b 2 ) + CO3, b 2 )(f) s Jx 2 , b 2 



r n r p (2 - x 2 - x 3 )(jf 1T (x 3 , 6 2 )0*(x 2 , b 2 ) + rvr p (2 + x 2 + x 3 )^(x 3 , b 2 )(p s Jx 2 , b 2 
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hj(xi, x 2 , x 3 , h, b 2 ) exp[-S ef (t 2 f )]]j 



M, 



up 



Mi 



(26) 
(27) 



For the factorizable annihilation diagrams (g) and (h) 



ap 



ap 



-F n 



(28) 
(29) 



If the p meson replaced by u meson in Figure 1 and 2, the formulas will be the same, except 
replacing f p by and (f> p replaced by cj) u . 

In the language of the above matrix elements for different diagrams eq.( |i"2"|]2"9| ,), the decay 
amplitude for B° — > n + p~ can be written as 



M(B C 



IT p 



ep 



M^i + C2)-6(C 4 



-c 3 + c 



10 



-a 



+ M ep [{ 
+ M a 



£uC 2 — it ( C*4 ~ Cq + -Cg + Cio 



C3 + C4 — C*6 — Cg — -Cg — 2^* 10 



c 5 - \c 7 



+ F a 
3 



Ci + ^C 2 



-3C3 - c 4 - -C 7 - -c 8 

-C5 + C*6 — -C7 — -Cg 

3 2 



:C fl + C 



10 



(30) 



where £ u = V* b V u d, it — V-foVtd- The C-s should be calculated at the appropriate scale t using 



equations in the appendices of ref. |L0|- The decay amplitude of the charge conjugate decay 
channel 5° — > p +r K~ is the same as eq.(|30|) except replacing the CKM matrix elements C, u to £* 
and it to £ t * under the definition of charge conjugation C\B°) = —\B°). 
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The decay amplitude for I? — >• it p can be written as 
- 2M(B° 



-°P°) 
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The decay amplitude for B + — > p + 7r° can be written as 
V2M{B + ^p+7r°) = (F e + 2F a 



6 ( ^C 3 + C 4 4 
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The decay amplitude 
V2M(B + ^tt + p°) 



for S + — > Tr + p° can be written as 
1 
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(34) 



From eg. fl30| - |3"4| ) , we can verify that the isospin relation 



M(B° -> 7r + p-) + M(B° -> 7r-p + ) - 2_M(5° -> ttV°) 



(35) 



\/2A^(5 + -> 7r°p + ) + V2M(B + -> 7r+p°), 



holds exactly in our calculations. 

The decay amplitude for S + — > 7t + c<j can also be written as expressions of the above Fi and 
Mi, but remember replacing / p by f w and p by </> u . 



V2M(B + ^n + u) = F e [e„ (Ci + ^C 2 ) 
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(36) 



The decay amplitude for B 



•> tt°u; can be written as 
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4 Numerical calculations and discussions of Results 



In the numerical calculations we use 



Agr 4 ) = 0.25GeV, f n = 130MeV, f B = 190MeV, 

m = lAGeV, f p = f u = 200Me^, fj = fl = 160MeV, 
M B = 5.2792GeK, M w = 80.41Ge^. 



(38) 



Note, for simplicity, we use the same value for f p (f p ) and f u And this also makes it 

easy for us to see the major difference for the two mesons in the B decays. In principal, the 
decay constants can be a little different. For the light meson wave function, we neglect the b 
dependence part, which is not important in numerical analysis. We use wave function for </>^ 
and the twist-3 wave functions (p^ and <f>\ from (l(| 

3 

u 



2V6 



2V6 



f n x(l -x)[l + 0.44C 2 3/2 (2x - 1) + 0.25C 4 3/2 (2x - 1] 
1 + 0.43C 2 1/2 (2x - 1) + 0.09C 4 1/2 (2x - 1) 
1 - 2x) \l + 0.55(10x 2 - lOx + 1) 



(39) 
(40) 
(41) 



The Gegenbauer polynomials are defined by 



C\ /2 {t) = i(3t 2 - 1), C\ r \t) = i(35t 4 - 30i 2 + 3) 



1/2/ 



C 3 2 /2 (t) = §(5t 2 - 1), Ct r \t) = f (21t 4 - 14t 2 + 1) 



t3/ 



15 , 



(42) 
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whose coefficients correspond to mo = 1.4 GeV. In B — > 7rp, tto; decays, it is the longitudinal 
polarization of the p and u meson contribute to the decay amplitude. Therefore we choose the 
wave function of p and u meson similar to the pion case in eq. (|39| , ^TD [fTTH 

3 



4> p (x) = (p^x) 
<p s Jx) = 4>' u (x) 



x 



0.18C 2 3/2 (2x 



^= {3(2x - l) 2 + 0.3(2x - l) 2 [5(2x - ] 
+0.21[3 - 30(2a; - l) 2 + 35(2a; - l) 4 ]} 



— f T (l 



2x1 



1 + 0.76(10a; 2 - lOx + 1) 



(43) 

(44) 
(45) 



Here again, for simplicity, we use the same expression for p and ui mesons 
For B meson, the wave function is chosen as 

iV B „ ^ .„ f M 2 B x 2 



2V6 



f B x (1 - x) exp 



2ul 



(46) 



with u>b = 0.4 GeV. N B = 2365.57 is a normalization factor. We include full expression of 
twist-3 wave functions for light mesons unlike B — » titt decays |10[ . The twist-3 wave functions 
are also adopted from QCD sum rule calculations [T6| . These changes make the B — > p form 



factor a little larger than the B — >• 7r form factor [[TBI . However, this set of parameters does not 
change the B° 7r + 7r~ branching ratios. And we will see later that, this set of parameters will 
give good results of B — ► np and ttu decays. With the above chosen wave functions, we get the 



corresponding form factors at zero momentum transfer from eq. (|l3l , |2"3"D : 



0.30, 



A*^p = 0.37. 



They are close to the light cone QCD sum rule results [13]. 



The CKM parameters we used here are [19 



\V ud \ = 0.9735 ± 0.0008 |Kb/Kh| = 0.090 ± 0.025 
\Vcb\ = 0.0405 ± 0.0019 \V t * b V td \ = 0.0083 ± 0.0016. 

We leave the CKM angle 02 as a free parameter. 2 's definition is 



(47) 



02 = arg 



v ud v: b 



(48) 
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In this parameterization, the decay amplitude of B — > up (or hd) can be rewritten as 



M 



V* b V ud T - VttVuP 



v: b v ud T 



tb vtd* 



ze 



i{<f>2+S) 



(49) 



where z 



v* b v ud 



and 5 is the relative strong phase between tree (T) diagrams and penguin 
diagrams (P). z and S can be calculated from PQCD. The corresponding charge conjugate decay 
mode is then 



M 



v ub v: d T - v tb v; d p 



1 



ze 



= v ub v: d T 

Therefore the averaged branching ratio for B — > irp is 

Br = (\M\ 2 + |AT| 2 )/2 



i(-<h+S) 



(50) 



V ub V* d T | 2 1 + 2z cos 02 cos 5 + z A 



(51) 



where z 







p 


Vu b v: d 




T 



Eq.(^TJ) shows that the averaged branching ratio is a function of cos 02- 
This gives potential method to determine the CKM angle 02 by measuring only the averaged 
non-leptonic decay branching ratios. In our PQCD approach, the branching ratios and also the 
other quantities in eq. (|51|) are all calculable, such that cos 02 is measurable. However, there 
are still uncertainties in the input parameters of our approach as discussions below. More 
experimental data from BABAR and Belle can restrict these parameters in the near future. 

More complicated, there are four decay channels of B°/B° — > 7r + p~, B°/B° — > p + n~ . Due 
to BB mixing, it is very difficult to distinguish B° from B°. But it is very easy to identify 
the final states. Therefore we sum up B°/B° — > ir + p~ as one channel, and B°/B° — > p + 7i~ 
as another, although the summed up channels are not charge conjugate states. We show the 
branching ratio of 5°/^° -> n+p-, B°/B° -> p+n', B + -> n + p°, B + -> p+vr , and B + -> n + u 
decays as a function of 02 in Figure 0. The branching ratio of B°/B° — > ir + p~ is a little larger 
than that of B°/B° — > 7r~p + decays. Each of them is a sum of two decay channels. They are 
all getting larger when 02 is larger. The average of the two is in agreement with the recently 



measured branching ratios by CLEO and BABAR [2C] 



Bi(B 



o 



Br(B° -> n + p- + n-p + ) = 27.6t 8 7 i ± 4.2 x 10" 6 , 

= 28.9±5.4±4.3 x 10" 6 , 
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Figure 3: Branching ratios (10~ 6 ) of B°/B° — > n + p~ (solid line), B°/B° — ► p + 7r~ (dotted line), 
B + — > 7r + p° (dashed line), 5 + — > p + 7r° (dash-dotted line), and _B + — > 7r + Ci; (dash-dotted-dotted 
line), as a function of CKM angle 02- 



There are still large uncertainties in the experimental results. Therefore it is still early to fully 
determine the input parameters and to tell the CKM angle 02 from experiments. 

The most uncertain parameters in our approach are from the meson wave functions. In 
principal, they can be only restricted by experiments, namely, semi-leptonic and non-leptonic 
decays of B mesons. Our parameters chosen for the numerical calculations in eq. (|38|j46|) are best 
fit values from B — > -kit decays [ID], B — > irK B — > it, B — > p semi-leptonic decays and 
some other experiments. As in these decays, we show the uuj, dependence of Branching ratios 
Br(B°/B° — > 7r ± p =F ) in Figure |j. The dotted, dashed and dash-dotted lines are for uo^ = 0.36, 
0.40 and 0.44, respectively. They are also shown as a function of CKM angle 02- The two 
horizontal lines in the figure are BABAR measurements of la. From the figure, we can see that 
the branching ratio is quite sensitive to the Ub parameter. Fortunately, this parameter is also 
restricted from semi-leptonic decays J03]. In the near future, it will not be a big problem for 
us. 

In Figure |5|, we show the branching ratio of B°/B° — > 7r ± p :F decays: mo = 1.3 GeV (dotted 
line), m = 1.4 GeV (dashed line) and tuq = 1.5 GeV (dash-dotted line) as a function of 
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Figure 4: Branching ratios of B°/B° — ► 7r ± p :F decays: u;;, = 0.36 (dotted line), u>b = 0.40 
(dashed line) and u>b = 0.44 (dash-dotted line) as a function of CKM angle <p2- The two 
horizontal lines are BABAR measurements. 

CKM angle 02- m o is a parameter characterize the relative size of twist 3 contribution from 
twist 2 contribution. It originates from the Chiral perturbation theory and have a value near 1 
GeV. Because of the Chiral enhancement of mo, the twist 3 contribution is at the same order 



magnitude as the twist 2 contribution. Thus the branching ratios of Br(B°/B° 



7T ± p T ) are 



also sensitive to this parameter, but not as strong as the cj& dependence. 

The branching ratios of B + — > 7r + p° and B + — > 7i + u> have little dependence on 2 . They are 
a little smaller than the CLEO measurement |IJ showed below, but still within experimental 
error-bars. 



Br(B + -> ttV) = 10.4+3;2 ± 2.1 x 10" b , (54) 
Br(B + -> ti + lu) = 11.3±|;| ± 1.4 x 10" 6 . (55) 
However, the recent BABAR measurement is in good agreement with our prediction for B + — > 

Br(B + -> tt+cu) = 6.6t?i ± 0.7 x 10" 6 , (56) 

where the error-bars are also smaller. The preliminary result of Belle shows that the branching 
r+ p° is around 6 x 10~ 6 p2| . This agrees with our prediction in Figure |[ 
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Figure 5: Branching ratios of B°/B° — > 7r ± p =F decays: mo = 1.3 GeV (dotted line), mo = 1.4 
GeV (dashed line) and m = 1.5 GeV (dash-dotted line) as a function of CKM angle 02- The 
two horizontal lines are BABAR measurements. 

The averaged branching ratios of B° — > 7T°p° and tc°uj are shown in Fig.[]. They also have a 
large dependence on 02. The behavior of them is quite different, due to the different isospin of 
p° and uj. But their branching ratios are rather small around 10~ 7 . They may not be measured 
in the current running B factories, but may be possible in the future experiments, like LHC-B 
and NLC. The recent BABAR upper limit of the channel is |2(^ 



Br(B° -> 7r°p°) < 10.6 x 1(T 6 . 

This is still consistent with our predictions in standard model (SM). 
Using eq. (|49| - pOD , we can derive the direct CP violating parameter as 

\M\ 2 - \M\ 2 



(57) 



a dir 
A CP 



\M\ 2 + \M\ 2 

2 sin <p2 sin 5 



(58) 



1 + 2z COS 02 cos 5 + z 2 

Unsurprisingly, it is a function of cos 02 and sin 02. They are calculable in our PQCD approach. 
The direct CP violation parameters as a function of 02 are shown in figure [7[ The direct CP 
violation parameter of B + — >• 7r + p° and B° —>■ 7r°p° are positive and very large. While the 
direct CP violation parameter of B + — ► p + ir° and B° — ► tt lj are negative and very large. The 
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Figure 6: Branching ratios (10 7 ) of B° — > 7r°p° (solid line), B° — »■ n°Ld (dotted line), as a 
function of CKM angle <p2- 

large strong phase required by the large direct CP asymmetry is from the non-factorizable and 
annihilation type diagrams, especially the annihilation diagrams. This is the different situation 



in Factorization approach where the main contribution comes from BSS mechanism [23fl and 



the annihilation diagram has been neglected The direct CP violation of B + — > tx + uj is 



rather small, since the annihilation diagram contributions in this decay is almost canceled out 
in eq . (p6|) . The preliminary measurement of CLEO shows a large CP asymmetry for this decay 
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A CP {B + -> 7i + u) = 34 ± 25%. 



(59) 



Although the sign of CP is in agreement with our prediction, the central value is too large. If the 
result of central value remains in future experiments, we may expect new physics contributions. 

For the neutral B° decays, there is more complication from the B° B° mixing. The CP 
asymmetry is time dependent : 



A CP {t) ~ Aq P cos(Amt) + a e+e i sin(Amt), 



(60) 



where Am is the mass difference of the two mass eigenstates of neutral B meson. The direct 
CP violation parameter A^ p has already been defined in eq.(^). While the mixing-related CP 
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Figure 7: Direct CP violation parameters of B + — > h + oj (solid line), B + — > 7r + p° (dotted line), 
B + -> p + 7T° (dashed line), 5° -> p ir° (dash-dotted line), B° -> tt°uj (dash-dotted-dotted line), 
clS db function of CKM angle 02- 



violation parameter is defined as 

-2Im(X CP ) 



^ 1 + \Xcp 

where 



2 



(61) 



V t tV td (f\H eff \B") 



\ _ " tb " ta\J l^ejji^ / 



Using equations (E&EO), we can derive as 



i I 7p i{8-<f> 2 ) 

\cp = , ^ v (63) 

Acp and a e + e ' are functions of CKM angle 02 only. Therefore, the CP asymmetry of B — > 7rp 
and 7ro; decays can measure the CKM angle 02, even if for the neutral B decays including the 
BB mixing effect. 

If we integrate the time variable t, we will get the total CP asymmetry as 

Acp = — !-5^& + T^W. ( 64 ) 
1 + x z 1 + x A 



with x = Am/r ~ 0.723 for the B° — B mixing in SM [fl9 1 . The total CP asymmetries of 
B° — > 7r°p°, 7r°u; are shown in Figure ^. Although the CP asymmetries are large, but it is still 
difficult to measure for experiments, since their branching ratios are small around 10~ 7 . 
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Figure 8: Total CP asymmetries of B° — > 7r°p° (solid line), and B° — > iy°u (dotted line), as a 
function of CKM angle <f) 2 - 

The CP asymmetries of B°/B° — > 7r ± p =F are very complicated. Here one studies the four 
time-dependent decay widths for B°(t) — * ir + p~ , B°(t) — > 7T~p + , B°(t) — > 7r~p + and -B°(t) — > 



7r + p |2g, |23, [28|]. These time-dependent widths can be expressed by four basic matrix elements 

g = (n+p-\H eff \B°), h = («+p-\H eff \B°), 

(65) 

g = (n-p+\H eff \B°), h = (7r'p+\H eff \B°), 

which determine the decay matrix elements of B° — > ir + p~ & ir~p + and of B° — > n~ p + & Tc + p~ 
at t = 0. The matrix elements g and h are given in eq. ( |30| , |3T| ) . The matrix elements h and g 
are obtained from h and g by changing the signs of the weak phases contained in the products 
of the CKM matrix elements. We also need to know the CP-violating parameter coming from 
the B° - B° mixing. Defining: 

B v = p\B°) + q\B°), 

B 2 = p\B°)-q\B°), (66) 



with \p\ 2 + \q\ 2 = 1 and q/p = y H 2 i/ Hi 2 , with = — i/2Tij representing the \ AB\ = 2 
and AQ = Hamiltonian. For the decays of B° and -B , we use, 



P v tb v t * d 
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e-** 1 . (67) 



So, \q/p\ = 1, and this ratio has only a phase given by — 20i. Then, the four time-dependent 
widths are given by the following formulae (we follow the notation of 



r(B°(t) - 


-> 7T~ 


V) 


= e" 


-id* 


2 + |h| 2 ){l + 


a e / 


cos Ami + 




_ e / sin Amt} , 


T(B°(t) - 


-> 7T 


> + ) 


= e" 




2 + |fc| 2 ){l- 


a f / 


cos Amt — 




sin Amt} , 


T(B°(t) - 


-> 7T 


> + ) 


= e 




2 + |ft| 2 ){l + 


<V 


cos Amt + 
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(68) 



where 



a,/ 



10 


2 _ 


h\ 2 


1.9 


2 + |/l| 2 



\p 9 



e+e i + IW 



(69) 



|/i| 2 + |#| 2 " 1 + |^| 2 
We calculate the above four CP violation parameters related to B°/B° — > 7r ± p T decays in 

PQCD. The results are shown in Fig.^J] as a function of 02- Comparing the results with the 

factorization approach p4| , we found that our predicted size of <v and a? are smaller while 

a e+e i and a e+ g' are larger. By measuring the time- dependent spectrum of the decay rates of B° 

and B°, one can find the coefficients of the two functions cos Amt and sin Amt in eq.(|68|) and 

extract the quantities <v, a e+e ', a?, and a e+ ?. Using these experimental results, we can tell the 

size of CKM angle 02 from Figfi 



5 Summary 

We calculate the B° -> 7r+p", B° -> p+vr^, 5+ -> p+7r°, 5+ -> 7r+p°, 5° -> 7r°p°, 5+ -> 
7r + u; and B° — > tt uj decays, together with their charge conjugate modes, in a perturbative 
QCD approach. We calculate the B — > 7r and i? — > p form factors, which are in agreement 
with the QCD sum rule calculations. In addition to the usual factorization contributions, we 
also calculate the non-factorizable and annihilation diagrams. Although they are sub-leading 
contributions in the branching ratios of these decays, they are not negligible. Furthermore 
these diagrams provide the necessary strong phases required by the direct CP asymmetry 
measurement. 
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Figure 9: CP violation parameters of B°/B° — > n ± p T decays: a e > (solid line), a? (dotted line), 
a t+e i (dashed line) and a e+ e' (dash-dotted line) as a function of CKM angle fa. 

Our calculation gives the right branching ratios, which agrees well with the CLEO and 
BABAR measurements. We also predict large direct CP asymmetries in B + — > p + 7r° and 
B + — > 7r + p° decays. Including the BB mixing effect, the CP asymmetries of B° — > tt°uj and 
B° — > 7T°p are very large, but their branching ratios are small in SM. The CP asymmetry 
parameters of B° — > 7i + p~, B° — > p + 7i~ decays require the time dependent measurement of 
branching ratios. 
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A Related functions defined in the text 

We show here the function h^s, coming from the Fourier transform of H^°\ 

h e (xi, x 2 , h, b 2 ) = K (y/x 1 x 2 m B b 1 ) [9{bi - b 2 )K {^ l /x^m B b l ) I (^x~ 2 m B b 2 ) 

+6{b 2 - h)K {y/x~ 2 m B b 2 ) J (v&i B &i)] S t (x 2 ) (70) 
h d (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = a s (t d )K (-iy/x 2 x 3 m B b 2 ) 

x [9(bi - b 2 )K ( v /xix 2 m jB 6i) I {^xix 2 m B b 2 ) 
+9(b 2 - bi)K (^xix 2 m B b 2 ) I (V^i^raB&i)] (71) 
tif(xi,x 2 ,x 3 ,b 1 ,b 2 ) = K (-iy/x 2 x 3 m B bi) 

x [9(bi - b 2 )K {-i^x 2 x 3 m B bi) J (y/x 2 x 3 m B b 2 ) 
+9(b 2 - b^Ko {-i^x 2 x 3 m B b 2 ) J (y/x 2 x 3 m B b 1 )] (72) 
hj(x 1 ,x 2 ,x 3 ,bi,b 2 ) = K (yx~ 2 + x 3 - x 2 x 3 m B b^j a s (t 2 f ) 

x [9(bi - b 2 )K (-i^x 2 x 3 m B bi) J (^x 2 x 3 m B b 2 ) 
+9(b 2 - h)K (-i^/x 2 x 3 m B b 2 ) J (y/x 2 x 3 m B b 1 )} (73) 
h a (x 1 ,x 2 ,b 1 ,b 2 ) = K {-iy/xxX 2 m B b 2 ) S t {xi) 

x [9(bi - b 2 )K (-iy/xim B bi) J {y/xlm B b 2 ) 

+ 9(b 2 - bi)K (-iyfx[m B b 2 ) J (y/x~[m B bi)} , (74) 



where J is the Bessel function and K , 1 are modified Bessel functions K (—ix) = — (n/2)Y (x)- 
i(ir /2)J (x). The threshold resummation form factor S t (xi) is adopted from ref.[|i"3"[ 

St{x) = 2 -^[T^ [x{1 - x)] ^ (75) 

where the parameter c = 0.3. This function is normalized to unity. 
The Sudakov factors used in the text are defined as 



S a b(t) = s (xim B /V2, bij + s (x 2 m B /V2, 6 2 ) + s ((1 - x 2 )m B /V2, k 
1 



Pi 



ln Inft/A) ln lnft/A) 



■ln(6iA) -ln(6 2 A) 
S c d(t) = s [xim B /V2, b^j + s (x 2 m B /V2, b 2 J + s ((1 - x 2 )m B /V2, k 

+s (x 3 m B / y/2, b^ + s ((1 - x 3 )m B / V2, b^j 



(76) 
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1 

Ti 



21nJm + h, ln ^ 



ln(6iA) ' *"-ln(6 2 A) 
S ef (t) = s (xxtub/VQ,, hj + s (x 2 m B /V2, 6 2 ) + s ((1 - x 2 )m B /V^, k 

+s (x 3 m B / y/2, 6 2 ) + s ((1 - x 3 )m B / y/2, 6 ; 
1 



(77) 



lnJ« + 21n. ln ^ 



ln(fciA) -ln(6 2 A) 
S g h(t) = s (x 2 m B /y/2, 61) + s (x 3 m B /V2, 6 2 ) + s ((1 -x 2 )m B /y/2, &j 

1 

"A 



+ s I-X3 m B /V2,6 ; 



ln ln(t/A) + ^ ln(t/A) 



(78) 



(79) 



ln(fciA) -ln(6aA)_ 

where the function s(q, b) are defined in the Appendix A of ref.pO[|. The scale tj's in the above 
equations are chosen as 



u 

,2 



max( v /i 2 "m B , l/b u l/b 2 ) , 

max(y/xlm B , 1/bi, l/b 2 ) , 

max( v /i 1 i 2 m B , y/x 2 x 3 m B , l/b u l/b 2 ) , 

max( v /x 2 x 3 m jB , l/b 2 ) , 

max(^x 2 x 3 m B , + £3 - x 2 x 3 m B , 1/bi, l/b 2 ] 



(80) 
(81) 
(82) 
(83) 
(84) 
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